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We discuss the evolution of cosmological tensor perturbations in the RSII model. In Gaussian 
normal coordinates the wave equation is non-separable, so we use the near-brane limit to perform 
the separation and study the evolution of perturbations. Massive excitations, which may also mix, 
C**) , decay outside the horizon which could lead to some novel cosmological signatures. 
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^ ■ I. INTRODUCTION 

The brane-world idea, according to which the observable universe is a hypersurface (brane) embedded in a higher- 
^v^j ■ dimensional spacetime (bulk), has attracted the attention of particle physicists and cosmologists in recent years — 
for reviews see 0, H H Q ■ The idea is motivated from developments in string theory and M-theory, although many 
t-H , of the models which have been developed are phenomenological. Two such models were constructed by Randall and 
Sundrum [1|||, hi which the bulk spacetime is a five-dimensional Anti-de Sitter (AdS) spacetime with a small length- 
scale (~ 1 mm or less) . In this article, we will be considering the second |(| of these models which has only one brane 
with a large, positive, bare tension to balance the curvature of the AdS bulk. The cosmological background solutions 
qq \ in the bulk spacetime were calculated shortly after the model was proposed. 

Cosmological perturbations around this bac kgro und, however, have resisted a solution in the face of a considerable 
research effort ESIUHHElHEm BJP El El El El El HI M El El El H3 (For a review on 
cosmological perturbations in brane- worlds see |3l]].) Although diverse formalisms for cosmological perturbations in 
brane models exist, the evolution of perturbations during the different cosmological epochs is not yet fully understood. 
The problem is that the bulk equations are partial differential equations which are subject to spacetime dependent 
boundary conditions on the brane. Therefore, the full problem can cannot be reduced to a four-dimensional effective 
\ problem. In addition, there is a large difference in scales: the AdS length-scale is of the order of 1 mm or less, whereas, 
the perturbations have to be evolved over many orders of magnitudes of conformal time, a ratio in scales of about 10 26 . 
Because of this, the wave equation for the evolution of the perturbations is difficult to solve numerically, at least by 
brute-force numerical integration: the spatial length scale of 1 mm means that we would have to use a corresponding 
time scale as the time-step, meaning the code would have to be run for the order of 10 26 time steps! 

In order to test the brane-world idea against cosmological data, it is extremely important to understand the evolution 
of perturbations. One of the most powerful discriminating tools between cosmological models is the power spectrum 
of the cosmic microwave background (CMB) anisotropics, which is, effectively, a snapshot of the perturbations in the 
universe when it was about one thousandth of its current size. At about this time in the history of the universe, 
two related events occur: atomic recombination and last scattering of photons. Before atomic recombination, the 
temperature of the universe is so high that protons and electrons form a plasma, but during recombination they come 
together to form atomic hydrogen. At approximately the same time, the mean free path of the photons due Thomson 
scattering in the plasma increases rapidly from being very small to very large (effectively, from zero to infinity) meaning 
that most of the photons arriving at the earth from the cosmological plasma will not have scattered subsequently. Once 
foreground effects due to our galaxy and external galaxy clusters have been removed, the electromagnetic radiation 
seen from earth is an image of this surface of last-scattering, or, more precisely, the intersection of it with our past 
light-conc. However, because of perturbations in the early universe, this will occur at slightly earlier or later times 
at different points of space, so the CMB radiation will be slightly brighter or dimmer at different points on the sky. 
The upshot of this is that, by measuring the temperature fluctuations in the CMB, we can deduce the perturbations 
in the metric on the surface of last scattering. 

The details of the perturbations in the metric at the time of last-scattering and their subsequent evolution depend 
on the theory of gravity. In brane-world models, Einstein's equation becomes modified, at least at high energies, and 
it is plausible that these effects can propagate into the regime relevant to the creation of CMB anisotropics. Clearly, 
the evolution of cosmological perturbations in this model is of significant interest. 
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The modifications in the Randall-Sundrum (RS) model are twofold: first, in Einstein's equation a term quadratic 
in the energy-momentum tensor appears. In the context of cosmology, this implies that at high energies the Hubble 
parameter on the brane is proportional the energy density of matter p and not proportional to ^fp, as in General 
Relativity. However, in the RS model one can show that this term is negligible at low energies, i.e., at energies much 
smaller than the brane tension. The other correction which appears in Einstein's field equation is the five-dimensional 
Weyl tensor projected onto the brane. This term encodes the influence of the bulk gravitational field on the dynamics 
of the brane. In the case of a homogeneous and isotropic brane this term vanishes in the case of an AdS bulk, but is 
non-zero if the bulk contains a black hole (Schwarzschild-AdS). From the point of view of 4D cosmology, this term 
behaves like a smooth energy component with equation of state P — p/3 and was therefore dubbed "dark radiation" 0. 
Its response to perturbations is somewhat different |l7j . 

The purpose of this article is to study the evolution of tensor perturbations in the RS brane-world. Although 
these perturbations are easier to understand than the observationally interesting case of scalar perturbations, we 
will see that their evolution is already very complicated. There are several coordinates system which one might use 
to tackle the problem. In the Gaussian normal (GN) coordinate system the brane is at rest but the bulk metric 
components have a complicated time-dependence: the brane boundary condition are easy to impose, but the solution 
of the corresponding wave equation is much more subtle since it is non-separable. Alternatively, one can formulate 
the problem in a coordinate system where the brane moves, but the bulk spacetime is manifestly static: one can then 
solve the wave equation in the bulk, but the boundary condition on the brane is difficult to impose. In this paper 
we will discuss solutions in the near-brane limit of the GN coordinate system. These solutions are the most direct 
an alog ue of the well studied mode functions of Minkowski space. Some of the material presented here first appeared 
in |32j and ■ The fact that the equation of motion is not separable is likely to lead to mode mixing and this will 
be discussed in future work. 

The paper is organized as follows: In Section^ we discuss the different scales involved and discuss different ways 
to attack the problem. In Section fllll we present the setup and derive the bulk wave equation and the brane boundary 
conditions. We also discuss the issue of initial conditions. In Section ITVl we present the near-brane solutions in the 
different cosmological epochs. In Section we discuss the problem of initial conditions. We present our conclusions 
in Section IVll 



II. OVERVIEW 

For simplicity, we will only consider the second RS model 6] which has a single brane of constant positive bare tension 
in an Anti-de Sitter (AdS) bulk which has a reflection symmetry between the two sides of the brane. Furthermore, we 
will restrict our attention to the case of a spatially flat universe, which is observationally favoured. The bulk Einstein 
equation for this model is then 

2 

= g A 3ai „ . (I) 
The Gauss relation gives the four-dimensional Ricci tensor as 

Rnu = ^ + KK^ - K/K vp - nV . (2) 
where the extrinsic curvature tensor of the brane is defined as 

Kfj. u = -g/V p n v , (3) 
with n v being the unit covector normal to the brane, and 

HV = C^nPrf , (4) 
is the contribution from the Weyl tensor of the bulk spacetime . The junction conditions [34| become 

K a p = ^{ 37 "a/3 + (T - \)g a p} , (5) 
whilst the definition of K a /3, combined with the bulk Einstein equations Q, yields 

£ n K a0 = W a/3 + - - K aj Kf , (6) 

£ n g a f3 = —ZK a p . (7) 

We will study perturbations by linearizing these equations about the cosmological background solution. For simplicity, 
we will focus on the case of a spatially flat universe, i.e., where the surfaces of isometry are flat planes, and consider 
only the case of a pure AdS bulk where W a p = for the background. 
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A. The variety of coordinates and gauges 

There are a variety of formalisms available for studying this problem. In this section, we explore some of these 
possibilities and discuss their various merits and drawbacks. There is no clear favourite, but some gauges are better 
for studying certain aspects of the problem. 

There are two different natural choices of coordinates in which it is convenient to express the background, in addition 
to the possibility of a coordinate-independent formalism. One natural choice is the Gaussian normal (GN) coordinate 
system. In these coordinates the metric for the background (non-perturbed) spacetime has the form 

ds 2 = -n(r, C) 2 dr 2 + o(t, C) 2 <% dxW + d( 2 . (8) 

The main advantage of this coordinate system is that the brane remains at a fixed value of one of the coordinates, 
meaning that imposing the boundary condition due to the perturbed junction conditions is simple. The main disad- 
vantage is that the metric components a and n have complicated functional forms and, as a consequence, the equations 
of motion are not separable. 

The other natural coordinate systems are ones where the static nature of the background is manifest. Another way 
of writing the line element for a flat cosmology in an AdS bulk is 

ds 2 = r) a p dx a dx 13 = —(^-dt 2 + % dx i dx j + dz 2 j , (9) 

where the spacetime is manifestly conformally flat. (Here, I is the AdS length-scale.) This coordinate system makes 
the five-dimensional linearized Einstein equations simple. The price paid for this simplicity is that the brane is no 
longer at a fixed coordinate value, but has locus given by z = I /a, making the boundary condition much more difficult 
to impose. 

In addition to the various coordinate systems for expressing the background, there are also various gauge choices 
for expressing the perturbations, as well as gauge- independent formalisms such as those used in [Tol Ull 12a . l26l l3fj | . 
The gauge we will use is the GN gauge, which arises by requiring that both the perturbed and unperturbed metrics 
have GN form, so that the perturbed line element has the form 

ds 2 = -n 2 (l + <t>) dr 2 + 2anh dr dx l + a 2 (5 tJ + /,, , j d.r'dr' + d( 2 . (10) 

When using a GN background, this is a very natural choice and has the advantage of making it trivial to extract 
the perturbation of the brane metric, g M „. Note that this does not completely fix the gauge since it it possible to 
impose additional requirement on the values taken by certain components of perturbation variable on the brane. For 
example, we could impose (j) — — on the brane, which would make the perturbation of the brane metric that 
of the synchronous gauge, a gauge often used in the treatment of perturbations in the standard, four-dimensional 
cosmology 35]. 

If we chose to work in the conformally Minkowski coordinates of 10 , we would naturally write the perturbed metric 

as 

ds 2 = ^ (rj a p + ha/s) dx a dx? , (11) 

which would allow us to make various gauge choices. For example, [2(j used the transverse-traceless (TT) gauge 
where h a p satisfies 

rf^h afi = , d a h a = . (12) 

Late in the history of the universe, the motion of the brane in this coordinate system is slow, so that the time 
coordinate, t, is almost the same as conformal time on the brane and so the TT gauge is approximately a synchronous 
gauge. This setup has the advantage that the resulting perturbed Einstein equations exactly soluble. However, in 
this gauge the position of the brane is not on the same locus as for the background but is displaced, an effect which 
has been dubbed "brane-bending" in the literature. 



B. Physical scales in the problem 

We will be considering the evolution of perturbations by taking the Fourier transform in the three spatial directions 
of the surfaces of isotropy and then evolving each Fourier mode. Thus, each mode has associated with it a length-scale 
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given by where k is the magnitude of the Fourier transform variable. Physically, we interpret k as the comoving 
momentum of the mode and fc _1 as the comoving wavelength. There is another length-scale, -ff -1 , given by the time 
dependent expansion rate of the universe, equivalent to the Hubble radius of the Universe. In the standard cosmology, 
these are the only two length-scales, giving rise to a single dimensionlcss parameter, kH^ 1 , which depends on time. 
In the language of cosmological perturbation theory, this parameter tells us whether the mode is outside or inside the 
Hubble horizon. 

For perturbations in the RS model, there is another length-scale, namely the AdS length-scale, I. Thus, we have 
three dimensionless parameters: kH^ 1 , IH and kl, two of which are independent. The first two of these depend 
on time: fc-ff -1 becomes larger at later times, whereas, IH is smaller at later times. Both will, of course, remain 
approximately constant during a period of inflation. As in the four-dimensional case, kH~ l gives the scale of the 
mode relative to the horizon size. Typically, the AdS length-scale, / will be less than 1 mm whereas the length-scale 
of the perturbations fc _1 will be greater than 1 Mpc, so kl will be an extremely small number. Except in the very 
early universe, IH will also be small since it is the Hubble parameter measured in time units corresponding to 1 mm. 



III. SETUP 



As already discussed, there are a variety of formalisms, each with its own merits. Here, we shall use a GN coordinate 
system and gauge, so that the perturbed line element takes the form 



d S ^ 



-n 2 (l + 4>) At 2 + 2anbi drdx* + a 2 (Si 



(13) 



The brane is located at £ = 0. Note that our choice does not completely fix the gauge as, for example, we can set 
4> and bi to be zero on the brane, thereby having the brane metric in synchronous gauge. As a shorthand, we will 
denote derivatives with respect to r and ( by dots and dashes respectively. Note that and bi are irrelevant for the 
tensor perturbations studied here. 

The observed matter, T a p, will comprise a background part into a background part 



a — 



-p 


' 








(14) 



and a perturbation 



5T a « = 



-Sp 


an 1 (p + p)(vj + bj) 


—a^ 1 n(p + p)v l 


SpS^-YSj 



(15) 



with £% = 0. This is similar to the formalism used in [2]j. It has the great advantage that the brane metric in 
synchronous gauge can easily be read-off and standard CMB computer code (such as CMBFAST |35j ) can be applied to 
solve the Boltzmann equations and determine the power spectrum of matter and CMB anisotropics. 

For the RS model, the bulk spacetime is pure AdS, which has a length-scale, I, related to the cosmological constant 
by 



The form of the functions a(r, £) and n(r, £) was found in j7J to be 



o(t,{) = o(t,0) 



ill 



i i T ,C) = n ( T ,0 = 



U(r) 
6A 

d(r,Qa(r,0) 
d(r,0) 



6A 



cosh(2C/0 



-U(rf 



6A 



sinh(2C/0 



1/2 



(16) 

(17) 
(18) 



where U(t) = A + p(r) is the total energy density supported on the brane, i.e., the bare tension, A, plus energy density 
of observed matter, p(r) 48]. The metric component n(r, £) is chosen so that r corresponds to conformal time on the 
brane. We will also tune the bare tension against the bulk cosmological constant so as to make the effective brane 
cosmological constant zero, for which we need to take A = 3/(4ttG1 2 ). With this we get 



o(r,C)=o(r,0) [e-^-^ S inh(C/0 



(19) 
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For convenience, we list the non-zero Christoffel symbols for this background, which are 

rO _n -i a t -o _ aa 

J- 00 - - , Oj ~a~ ij ~ ^2 d H ' ^ U J 

and the non-zero components of the Ricmann tensor, which are 

= na-na §i ^ = ^ g[i 

The bulk Einstein equations give the following relations between the metric components 

2 , n" a an a'n' 

- o A = 3— + 3— + 3 , 22 

3 n an z an a an 

2. n" a" , N 

-»A = — + 3— , 23 
3 n a 

z , a an a n a a a , 

--A = 2+ — + + — - 2 ^ + 2 — - (24) 

3 an an° an a a z n z a 

a! an' 

= - + . (25) 

a an 

The brane, at £ = has extrinsic curvature given in Gaussian normal coordinates by 

1 d _ 

= ~2dC 3a/3 ' ( ' 26 ' ) 

which has components 

K m = nn , K i = , K i3 = -aa'S.^ , (27) 

from which we can deduce that 

K = -U a - + ^\. (28) 
\ a n J 

Note that some authors use the opposite sign convention for the extrinsic curvature. 

A. Equations of motion and brane boundary conditions 

Since we have chosen GN gauge, where the perturbed metric is in GN form, it is natural to consider the perturbation 
of the Gauss relation © which is 

SR^u = + SKK^ + K8K pv + K^K v °5g pa - 2K ( /SK v)p - SW pv , (29) 

and perturbing © and Q gives 

^SK a0 = 2K (a PSK 0)p - KJK p a 8g p<7 - ^6g a(3 - 6W a p (30) 

^Sg a0 = 2SK aP (31) 

substituting these into H29[) and equating with the perturbations of the Ricci tensor calculated from the metric gives 
the evolution equations for the metric perturbations. We will evaluate these evolution equations in the special case 
of tensor perturbations. Perturbing the Israel conditions of (JSJ) gives 

5K m = -|n 2 { 2Sp + 35p+(2p + 3p)<t>}, (32) 

8K 0i = ^«n{3(p + p)v, + (2p + 3p - A)6,} , (33) 

SKij = ~^a 2 ^SpS i3 + (p + X)h i3 + 3S y | , (34) 

which will provide the boundary conditions we need at the brane. 



6 



brane 




FIG. 1: Specification of boundary data: on the left, at the brane and the horizon; on the right, at the brane and on some initial 
surface. The bold lines show where boundary date is specified. The dotted lines are for guidance and indicate surfaces with a 
Minkowski four-metric. 

B. Initial conditions and boundary conditions at the horizon 

Since we are aiming to evolve the equations of motion over the semi-infinite interval < £ < oo it would appear 
that a another boundary condition is not necessary. However, the coordinates used do not span the universal covering 
set of the AdS spacetime (see [3(| for more details on the universal covering set) so we could expect that gravitational 
waves could cross the AdS horizon. The choice of boundary condition on this horizon is somewhat arbitrary since we 
can have no knowledge of the other side; this being a philosophical shortcoming of the Randall-Sundrum model. A 
reasonable, and popular, choice is to have no incoming radiation at the horizon, as advocated in [37j . 

However, it is not possible to impose this boundary condition in addition to arbitrary initial conditions on some 
time-slice. So we can impose boundary data in two different ways: 

1. Brane boundary condition and a condition at horizon. 

2. Brane boundary condition and initial conditions at some time. 

The second of these is more in keeping with the orthodox approach to cosmological perturbation theory. However, it 
requires knowledge of what created in the first place. These are not as different as they might seem, we can gain a 
better understanding of how they are related, and of when it is necessary to specify boundary data at the horizon by 
considering Carter-Penrose conformal diagrams of the spacetime, see Fig. ^ It is apparent from the diagrams that 
specifying data on the horizon is a limiting case of specifying data on a timelike initial surface. 

To illustrate this further, it is useful to consider the solution of the ordinary, wave equation in one spatial dimension, 

(35) 

with a boundary condition defy/ dx — at x — 0. The spatial part of the mode solutions are then of the form cos(o;x), 
multiplied by a time dependent part of the same frequency, oj. If we then impose a no incoming radiation condition, 
the solution is completely determined, and is constant, so it is not possible to impose initial data in addition. This 
can be seen more clearly by writing the general solution in the form f(x — t) + g[x + t). The no radiation condition 
would force us to take g = and the brane boundary condition necessitates that / = constant. 

For studying the initial perturbation spectrum created in the inflationary era, it is useful to use the no-radiation 
boundary condition, as was done in ;3?j. Once we have the perturbations at the end of inflation, we can evolve them 
forward through the radiation eras, for which we will have initial data on a time-slice, and then through the matter 
era. 



C. Tensor modes 

The perturbations in the four-dimensional part of the metric naturally split into scalar, vector and tensor parts 
with respect to the Fourier transform variable, k % . To linear order in perturbation theory, they do not mix. Of these 
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the scalar perturbations are the most important since these contain the underdensities and overdensities giving rise to 
large-scale structure and so necessarily exist. As we have already suggested, it is much easier to solve the equations 
for the tensor perturbations. We have formulated the equations of motion for all three of these but will only provide 
solutions for the tensor modes. Some of the lessons learned should assist in tackling the more difficult problem of the 
scalar modes. 

The variation of the purely spatial part of the Ricci tensor is given by SRij = 8R°ioj + 5R k ikj, so the tensor part 
is given by 

-t, a 2 -rp a 2 ( a h\ ■ ™ a 2 I a ah a 2 \ T k 2 ^ , . 

lJ In 2 13 2n 2 \ a n ) lJ n 2 \a an a 2 ) J 2 y y ' 

where the superscript T represents either of the two tensor modes, which will evolve separately, to linear order, and 
will both obey the same equation of motion. Equations l|3U|) and give us 



*S = *£(°X). (37) 



1 d 

2dC 

= " -'X - - \^ (a 2 hl) . (38) 



Substituting these into the tensor part of (|29|l gives 



The purely spatial part of the Einstein equation Q evaluated for the background metric gives 



I , a „ a an a a n a 



-A = — r + 2^^r ? 2^- (40) 



3 an 2 a 2 n 2 an 3 a an a 
using this and equating and (|3*5|l gives 



a n \ ; ^ , r,n 2 , ^ „ / a' n' \ d , rp „ d 2 



^ + + ^ - \ 3 -a + i) 5r ft 5 - " 2 ^"S " < 41 > 

in agreement with the result in |18|. 

The boundary condition on hjj at the brane is given by the perturbation of the purely spatial part of the junction 
condition (|34(l which becomes 



= . (42) 



The tensor part can be written in terms of two polarization tensors ef, and e^- as 

f$ = h+e+ + h*e*, (43) 

so, in subsequent equations, we will simply write h to represent either h + or h x . 

For most of the history of the universe, there is effectively no matter source for the tensor modes so this condition 
will reduce to h! = on the brane. Only when the CMB photons develop a quadrupole moment in the late universe 
is this assumption no longer valid. 

The wave equation (|41(l is not separable and is difficult to solve in general, due to the complicated form of the 
scale factors a and n, see l|19[) . We will first study two cases where the equations can be solved exactly: the case 
of a Minkowski brane and the case of a de Sitter brane. The first of these was the case considered by Randall and 
Sundrum in their orig inal paper [f|; the second is very useful in cosmology for modeling an inflationary era in the 
early universe |18|, |37j • 



IV. BASIS SOLUTIONS 



A. Near-brane limit 



Our main interest is in cases which correspond to different cosmological eras. In an inflationary model, the universe 
will undergo a phase of accelerated expansion where the brane metric will be approximately de Sitter. After a 
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transition period, there will be an era where the matter content of the universe is dominated by a radiation fluid 
component, having equation of state p — p/3. There will then be another transition to evolution dominated by 
pressureless matter. Shortly after this transition, photons decouple from the cosmological plasma, producing the 
CMB radiation. Since the eventual aim of studying perturbations is to determine theoretical CMB power spectra for 
RS models and compare these to the observed CMB spectrum, we should try to solve the problem in each of these 
eras. 

Unfortunately, the equation of motion is not separable in the radiation and matter dominated eras so we will 
approximate the solution by one valid near the brane, ignoring the sinh(£/Z) term in (|19f) . The range of £ for which 
the near-brane limit is a good approximation depends on the density of the cosmological fluid. The value of £ for 
which a(r, £) = is given by 

< = M 1 + 7)' (44) 

which will give us a qualitative idea of the domain of validity of such a solution. Since p will decrease with time in the 
radiation and matter eras, this domain of validity will get wider later in the history of the universe. The near-brane 
limit can also be thought of in some sense as a low energy limit, but the fact that it is an approximation suggests 
that discussion of a low-energy effective theory needs to be careful thought out. 

Let us see where the locus of a = intersects the past light cone of a point on the brane today since this will 
correspond to the limit of causal validity of the near-brane limit. It is simpler to work in proper time on the brane, 
t, so that a null geodesic satisfies 

^=±n(t,C)^±e-^ 1 , (45) 



near the brane. This integrates to give 



e<' 1 = 1 + , (46) 



for the past light-cone, where to is the time today. We want to compare this to (|44|l to see when we can argue that 
the approximation is good on causality grounds alone. During a matter dominated era of the universe, the critical 
energy density is given by 

and the bare tension on the brane is related to the AdS length-scale by A = 3/(47rGZ 2 ). So, in the matter dominated 
era, (|44|l becomes 

/ t 2 \ 1/2 t 
e C/« = (l + 9 *.j « 3 * ; (48) 

which intersects the light-cone when t ss to/4- Of course, the approximation used in Q45[) is good near the brane but 
breaks down as the null geodesic approaches the line l)44|l so this is only and order of magnitude estimate. Nonetheless, 
this means that we cannot use a causality argument to justify the approximation throughout the whole history — we 
would like to have obtained a time before the end of inflation as the answer! It seems clear that some mode mixing 
will take place. 

Finally, we should point out that the near-brane approximation is also equivalent to an expansion in powers of 
(HI) 2 (or p/X). Since it involves ignoring the second term in (|19[) . one could do this either because smh(£/7) is small, 
or because p/X cx [HI) 2 is small. 

B. Minkowski brane 

For the RS case where the projected metric on the brane is Minkowski, the functions a and n have the simple, 
separable form 

a(T,Q=n(T,Q = e-V l . (49) 
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This simplifies giving 



Or, equivalently, from l|41|l 



lv k 2 . /A 1\ o. Id 2 



i h+ n Y h = 2^ + ±)a 2 h+^^ h ). (50) 



/i + fc 2 / l = e- 2 ^ . (51) 

This equation has a constant solution, corresponding to a metric perturbation proportional to e~ 2 ^, which is the 
same as the solution found by Randall and Sundrum|6|]. They interpreted this as a graviton zero- mode localised on 
the brane, thus having the appearance of a four-dimensional gravitational perturbation. 
To find a set of basis solutions, we perform a separation of variables, using the ansatz 

h(T,Q=il>(T)4>(Q. (52) 

The equation of motion Il- r >()[i then separates into 

4>m + (k 2 + m 2 ) ip m = , 

C-y^ + ™ 2 e 2C/ V m = 0, (53) 

where m 2 is a separation constant. Because m 2 has dimensions of mass, modes with m 2 ^ are referred to as massive 
modes, but it should be noted that m 2 can be negative. The solutions are 

ipm(T~) — ccos (^J to 2 + k 2 tJ + d sin y\J m 2 + k 2 rj , (54) 
MO = C + De 4 ^ 1 , (55) 
MO = e 2 ^ Z 2 (mle</ 1 ) , (56) 

where c, <i, C and D are constants of integration and Z v represents a linear combination of Bessel functions of order 
v . (The boundary conditions will determine which particular combination to choose.) Randall and Sundrum [(| also 
showed that, in the Newtonian limit, the contribution from the massive modes is sub-dominant, and so gravity has 
approximately four-dimensional behaviour. 

C. In a de Sitter era 

During an era of inflation in the primordial universe, the density of matter is almost constant, an so the spacetime 
is approximately a de Sitter universe. For this reason it is useful to consider a brane-world where the brane metric 
takes de Sitter form. For a de Sitter brane, the matter density of the background solution is a constant so scale factors 
a and n are separable functions. We shall work in conformal time, where they have the functional form 

a(r,0=n(T,Q = a (T)A(0, (57) 

where oq(t) = a(r, 0) and ^4(C) is given by 



* -(!+&).-<" + &.<". (58) 



The wave equation 1)41(1 becomes 



h + 2— h + k 2 h = A 2 h" + AA!Ati . (59) 
a 

We can find a set of basis solutions to this by performing the same separation of variables used in the RS case. Writing 
h = iP(t)4>(^), the wave equation splits into two parts: 

4> m + 2— ip m + (k 2 + m 2 al) ip m = 0, (60) 
A' m 2 

C+4— ^ + _0 m = O , (61) 
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where m 2 is a separation constant as before. 

In a de Sitter era p is a constant so, in conformal time, we have a(r, 0) = — (H- ln fr)~ l where H ln f is the energy scale 
for inflation, giving the evolution of the tensors as 

0, (62) 

(63) 

The solutions for tp are 

= t 3/2 (cT a + dr- a ) , (64) 

for k = 0, and 

i/> = r 3 / 2 (c,J a (fcr) + dY a (fcr)) , (65) 

for k ^ 0, where a = ^/9 — (Am 2 /H 2 nf )/2 and c and d are constants of integration. The equation for <f> reduces to a 
hypergeomctric equation, which we can be solved exactly, but the solution is very complicated and not particularly 
informative. Let us state briefly the m — solution, which is 

4>o = C + De- 4 ^ 1 . (66) 

This de Sitter case has been studied extensively in the literature. It was found that the Kaluza-Klein modes are not 
continuous, but there is a mass gap between the zero mode and the first massive mode. If one requires that a is pure 
imaginary then m > (3/2)H{ n {. Above that mass the spectrum is continuous. In |18| it was argued that the tensor 
perturbation left after an inflationary era would mainly be comprised of the zero- mode solution, h — -0o(r)0o(C); an d 
this result was made more quantitative in j37l | where Bogoliubov coefficients were calculated. Thus, the perturbations 
produced from a period of inflation will be the zero-mode, h = iI>o(t)4>o(t), to good approximation. This has been 
extended to the case of any conformally flat brane metric in [3^ . 



4>n 



-Iprr. 
T 

4 T 



: (fc 2 + 

o 

m 



1pm = 



D. In the radiation era 



Let us now consider the evolution of the tensor perturbations in the radiation-dominated epoch. As mentioned 
earlier, we cannot solve the equations exactly and so we use the near-brane approximation. As for the de Sitter era, 
we work in conformal time, where the the scale factors a and n have the form 

o(r,C) = n(T,0=ao(r)e-^. (67) 

We are assuming that the density is small compared to the bare tension of the brane in order for the near-brane 
solution to be a good approximation, so we can approximate the Friedmann equation by the usual, four-dimensional 
one 

da\ 2 8vrG 4 

Tr) =— Pa ' (68) 

which, if we choose the origin of time appropriately, has solution a = At for some constant A. 

Our near-brane approximation allows us to separate variables, so we try to find solutions of the form h = ip(r)(f)(Q 
as before. The wave equation then splits into two equations of motion 

i'm + -V>m + [k 2 + m 2 A 2 T 2 ] ip m = 0, (69) 

T 

C-J<4 + ™ 2 e 2C/ Vm=0. (70) 
The solution of Ij70|l is independent of k and is given by 

MO = C + De«' 1 , (71) 
MO = eK ' 1 3>(m/e</') , (72) 
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where C and D are integration constants, and, as before, Z v represents a linear combination of Bessel functions of 
order v with the boundary condition determining the particular combination. Assuming we have no source for the 
tensor perturbations, the boundary condition is simply cfi' — 0, which gives us basis solutions 

MO = e 2C/i jyi(mOJ 2 (mZe^) - Ji{ml)Y 2 (mle^ 1 ^ j , (73) 

with J v and Y v being, respectively, Bessel functions of the first and second kind. When k = 0, the solutions to Il69() 
are 

Mr) = c+-, (74) 

T 

^m{r) = r- 1 ' 2 Z 1/A (\mAr 2 \ , (75) 



and when k ^ the solution for to = is 



, , , sin for cosfcr ,_„. 
Mr)=c + d , (76) 



where c and d are constants of integration. 

For to ^ 0, the solutions can be written in terms of parabolic cylinder functions as follows. Writing i/j = y/r and 
changing the independent variable to t = fcr, (|69[1 becomes 

d 2 y ( m 2 A 2 ^\ 



dt 2 \ k 4 

which is Weber's equation and has parabolic cylinder functions as solutions |42|, |4.. 

k 2 I 



±(l + i)yMAk-Hj , where v m = ^— - - - . (78) 
We can use the asymptotic form 

D„(z) -e- z2 / 4 z^, for |z|>l, |argz| < 3tt/4, (79) 
to approximate the solutions by 

e k 2 ir/(8mA) f ■ ik 2 ) 

Vm ~ (W)VV? 6XP l-2 mAr2 " W l0g 7 ' as ( 8 °) 

which is used in Appendix A to derive the orthogonality relation. If k = this can seen to have the same form as 
H74J) as r — > oo. The general solutions can then be written as a superposition of these modes of the form 



h = C + + C - +J2 D Vm (±(1 + 8 )V^r) {Yi(TO/)J 2 (TO?e^) - ^(mO^^k^') } 

± 



dm . 



(81) 

The evolution of the zero-mode is exactly the same as for the usual, four-dimensional cosmology based on General 
Relativity. For all values of k, the very heavy modes decay outside the horizon. 

Since special functions of this nature are not always easy to visualize, let us try to approximate to solution to l|77|l 
using the WKB method. First define the dimensionless parameter 

k 2 

6 = —A ( 82 ) 

mA 

The WKB approximation for e small recovers the asymptotic expansion given above. When e is large, we rescale the 
time coordinate by T — t/e — mAr/k, so that (|77|l is rewritten as 

1 ^ 1 (1 + T 2 )y = (83) 



dT 2 
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We will use the formulation of the WKB approximation as described in [39j | , where a small parameter S is introduced 
and an asymptotic solution of the form 



^v\- s jrs n s n (T)\ 



(84) 



is postulated and the S n are functions to be determined. The parameter 5 will be determined in terms of e by a 
distinguished limit, i.e., the choice for which a non-trivial answer is obtained. Substituting (|84|l into the differential 
equation gives 



(85) 



The term proportional to (e<5) 2 is the leading term on the left and so must match the right hand side. Hence, we see 
that S — 1/e is the distinguished limit and 



S = ±i\/l + T 2 , (86) 
The coefficients of the other powers of e give differential equations for the other S n : for the next three, we have 

3T 2 - 2 



Si 



T 



2(1 + T 2 ) ' 
The solutions to these equations are 

So 
Si 

s 2 
s 3 



So = ±i- 



(l + T 2 ) 5 /2 ' 



3T(2T 2 - 3) 
8(1 + T 2 ) 4 



(87) 



± % - (T\/i + T 2 + arcsinh(T)) , 



log(l + T 2 ) , 



=F 



i T(T 2 + 6) 
24(1+T 2 ) 3 / 2 
2-3T 2 



16(1 + T 2 ) 3 



(89) 
(90) 

(91) 



which will give an approximation to y m . Note that constants of integration have been omitted because they merely 
correspond to multiplying the solution by a constant factor. 

It is useful to consider the early and late time behaviour of this approximation. When T is large, the asymptotic 
expansion (|80|l is recovered. When T is small, the S n are approximately 



so the solution is 



So ~ ±iT , 



Vn 



Si~-\t 2 



4 16 



exp 



T 2 



1 



4e 2 



exp 



±iT e 



4c 



which is a good approximation for kr <^ e. We see that the solution which is bounded at the origin is 

1 



Vv, 



T 



exp 



r t 2 ( 












sin 






4e 2 A 







(92) 



(93) 



(94) 



which is flat near the origin and then begins to decay. 



E. In the matter era 



We now move on to the final phase of the evolution of cosmological perturbations when the universe is dominated 
by a pressureless fluid. As with the radiation era, we approximate the solution near the brane, where the scale 
factors a and n have the form (|67|l . Again, we are quite late in the history of the universe, so l|68|) is valid to good 
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approximation and we have a = At 2 if the origin of time is chosen appropriately. Making the ansatz h = i/}(r)<p(^), 
the equations of motion separate once again and the equation for ip is 

$ m + -1pm + [k 2 + m 2 A 2 T i ] <ip m = . (95) 

The equation for (j> and its solutions are the same as in the radiation era (fTO - )) . For k = 0, we can solve l|9"5|l exactly, 
giving 

Mr) = c+4, (96) 



For k ^ 0, the zero-mode can be found exactly 

ipo(r) — {kr sin(fcr) + cos(fcr)^ H — -^/c-rcos(fcT) — sin(fcr)^ . (98) 

We will study the modes with k 7^ and m ^ by constructing approximate solutions, which will be appropriate 
in various limits. First, we will define the dimcnsionlcss variable 

and two rescaled time variables t = kr and T = E~ x / 2 t = {mA/k) l / 2 T which we will use in the limits e< 1 and e> 1 
respectively. The equations of motion for ip are then 

and 

+ + (101) 



Case e < 1: 

Let us consider first the case where e is small and use ()100f> . We can then find a solution using the WKB 
approximation as discussed in the previous section. We try to find an asymptotic solution of the form l|84|) . which, 
when substituted into (|100|) gives 



£ p (E 2 + j E ^™ + 7T E ^ n + e 2 = -i 4 . (102) 

We see that the distinguished limit is when £ — S and that the S n are given by the differential equations 

Sl = ~t 4 , (103) 

2S S! + Sq + jSo = 0, (104) 

2S S 2 + S 2 + S 1 + jS 1 + l = 0, (105) 

71 — 1 . 

2S S n + SjSn-j + Sn-i + -S„-i = , (106) 
for the first four terms. Solving for these, we find 

s = ±lit 3 , ^ = io g (i- 3 ), ^ = t^, s 3 = -|r- ( 10? ) 
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Substituting the first three of these into (|84|l along with the expressions for e and t, the solution with the first 
correction term is 

^ = C 5/(mW) ( _ _&\ d e - 5 /( m ^ 2 ^) gin ( 1 mAr 3 _ (108) 

Y r 3 V 3 2mAr J T V 3 2mAr J K ' 

This solution is good for small k and is progressively better for larger values of r. We see that in the limit k — > this 
solution reverts to the form (|97|l as one would expect. 



Case £> 1: 



Now consider the opposite limit where £ is large using l|101|) which is in a form suitable for using the WKB method 
as formulated above. We get S = e -1 / 2 and the following differential equations, 



So = ±tVl + T\ 
2S* S'i+S'o+4T- 1 S'o=0, 
2S* S 2 + Sf + Si + 4T- 1 S* 1 = , 
2S S 3 + 2S\S 2 + S 2 + 4T~ 1 S 2 = , 

for the first four terms; the last three equations simplify to 



Si 

These have solutions 

±s = i 



T 3 



1 + T 4 T 



±S 2 



?T 4 



2iT 2 (l + T 4 ) 



5/2 



35T 8 + ?T 4 ^ 

T 3 (l + T 4 ) 4 



Si 
±S 2 

s 3 



T\fl + T 4 + ^e"/ 4 F (W 4 T, ij , 
log (t 2 (1 + T 4 ) 1/4 ^, 
iT 3 



5iT 3 



7T 4 + 2 



Z e — /4 



F 



(e i7r/i T, ij-E (e i7r/4 T, i 



2T 2 (1 + T 4 ) 3 



where 



F (e"/ 4 T, *) = y 



da; 



E 



1 + ix 2 
1 — ix 2 



d;i" 



(109) 
(110) 
(111) 
(112) 



(113) 

(114) 
(115) 

(116) 

(117) 

(118) 



are, respectively, elliptic integrals of the first and second kinds. These expressions are quite complicated, so let us 
consider the two limits where T is small or large. 

When T is small, the differential equations I|109I112[) are 



So sa ±i . 

so that the solution is approximately 



Si 



2 
T 



S 2 ^T 



V>(T) 



= 1/( E T 2 ) 



rp 2 



■ exp 



rp 2 ' 



oc 



■ exp 



2 

~J<3 ' 



±i I kr + — 
kr 



(119) 



(120) 



which should be a good approximation if 1 -C kr -C y 7 ?. This approximate solution for small T does not depend 
on m, which is a consequence of the fact that the term involving m is the equation of motion has a T dependence. 
When T is large equations (|109I112[) become 



So ~ ±iT 2 , 



Si 



3 

T ' 



S 2 ~T 



7z 

2T» ' 



S 3 



35 



(121) 



15 



giving us 



e 7/(2 £ T 10 ) 

ip(T) « -3 exp 

valid for kr 3> y^. 



i ex P' 



7fc 2 



\2m 4 A 4 T 10 



exp 



mAr 3 



2m 3 A*T 7 



(122) 



V. IMPOSING THE INITIAL CONDITIONS 



We are considering tensor perturbations here and since the matter on the brane has no tensor perturbation early in 
the universe, there is no matter source for these perturbations in the metric. So the boundary condition at the brane 
is h' — 0. The C dependence of the solution near the brane is the same in all eras, and its evolution is given by J7UJ), 
with solutions given in (|7T)> and i|72|). The derivative of </> m (C) is 



so, if the junction condition is imposed, the massive modes are 



m (£) = e 2C/l {Y 1 {ml)J 2 {mle c - /l ) - Ji(77iZ)Y" 2 (mZe c/z )) 



for which 

q>' m (Q = me 3C/l (Vi(mZ) J x (mle c/l ) 
The zero-mode is simple a constant 

M0 = c, 



Ji(ml)Yi(mle c/l ) 



(123) 



(124) 



(125) 



(126) 



when the boundary condition is imposed. 

Let us now try to impose initial conditions at some initial time, Tj, which will, in practice, be the transition between 
two eras. This initial data must satisfy the boundary condition and will be expressible as a superposition of the mode 
solutions 4> m (C) ■ This superposition is most easily calculated by considering the C derivative of the initial data. So 
let us write 



h'{nX)= A(m)<p' m (() dm = mA{m)e 3C ' l Cx{m-M /l )&rn, 
Jo Jo 



where C v {m]z) is the cylinder function 



C v {m; z) = Yi(ml)J v {mz) — Ji(rnl)Y v {mz) . 



(127) 



(128) 



We can deduce the A(m) from the orthogonality relation given in the appendix. 

The contribution from the zero-mode is just a constant, and this is, of course, the bit that is undetermined 
by considering h' . So we have a prescription for imposing any given initial conditions. The prescription given is 
somewhat academic, at least for the application under consideration here, because the spatial parts of the mode 
functions, (j) m (^), are the same in all eras, so the spectrum of modes present at the beginning of an era will be the 
same as the spectrum at the end of the previous era. In particular if inflation produces a perturbation which is well 
approximated by the zero-mode, as indeed it does according to |l8t l37|. then only the zero-mode is present during 
subsequent eras too. Thus the evolution of tensor perturbations is the same as in the standard, four-dimensional 
case, and hence the contribution towards CMB anisotropics is the same. The absence of sources is responsible for this 
similarity, tensor matter perturbations on the brane would source the massive modes and give an answer different 
from the standard cosmology. The analysis presented here should extend to the scalar perturbations, although there 
are some additional technical difficulties. The scalar perturbations are not only the dominant contributions to the 
CMB anisotropics but are sourced by scalar matter perturbations so there will very probably be some signature of 
the RS model detectable in the CMB power spectrum. 
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VI. DISCUSSION AND CONCLUSIONS 



In this article we have discussed some aspects of tensor perturbations in the Randall-Sundrum brane-world. The 
problem can be formulated in several coordinate systems, each of which has its shortcomings. In GN coordinates, 
the brane is static and the boundary condition can easily be imposed. However, the equations of motion for tensor 
perturbations are not separable. By formulating the problem in static bulk coordinates in which the brane moves, the 
bulk equation can be solved but the boundary conditions are difficult to impose on the (moving) brane. Therefore, 
an approximation scheme is necessary to make progress in both systems. We have presented solutions in the de Sitter 
era, the radiation and matter dominated era using the near brane approximation first discussed in [33, H3 • 

The near-brane approximation is akin to an adiabatic treatment of the motion of the brane in that the warp 
factor a(r, £) has the same £ dependence as in the standard Randall-Sundrum picture for a Minkowski brane, but is 
multiplied by the cosmological scale factor. In this lowest-order approximation, independent modes persist throughout 
the radiation and matter eras. If this approximation is universally valid, which we have argued against on causality 
grounds, and the initial spectrum includes only the zero mode created during inflation then the observed tensor 
spectrum will be exactly that of the standard cosmological model. 

The zero mode is equivalent to the standard case of General Relativity, but the massive modes have some novel 
features. In particular they can be seen to decay and oscillate outside the horizon, whereas standard massless 
gravitons only begin to oscillate at horizon crossing. One important consequence of this is that initial power spectrum 
of fluctuations in the massive modes cannot be scale invariant as k — > 0. There must exist a cut-off in the spectrum 
at very small values of k. If one thinks of the perturbations as being particles created quantum mechanically during 
inflation then one can create those with k <C m just on energetic grounds. Moreover, if massive modes exist and they 
decay outside the horizon, the photon quadrupole will increase without the need for Thomson scattering since it will 
be sourced by the decay of the super-horizon gravitational potential. This could lead to a significant polarization 
signature in the angular power spectrum on large angular scales. 

It is worth noting that the equations of motion for the massive modes which we have solved here are equivalent to 
those of ordinary massless gravity with a mass term added to the Lagrangian in an ad hoc way via an addition term 
of the form 



where m is the graviton mass and is the 4D metric perturbation. Although such a theory is believed to suffer 
from the van Dam-Veltman discontinuity |44j in the graviton propagator in going from m = to m > 0, this theorem 
will not apply in the RSII case where there is likely to be a spectrum of modes including the zero mode. 

It should be possible to incorporate non-adiabatic effects, by taking into account the non-separability of the equa- 
tions. We anticipate that this will lead to mixing of the near-brane modes discussed in this paper, as considered 
in ^(J . This is likely to lead to non-trivial time evolution of the total perturbation amplitude and interesting 
cosmological signatures. The solutions found in this paper can be used to gain insights into the full problem. The 
non-linearity in the wave equation induced next-order correction to the expressions (|17|l and (|18|) for the scale factors 
a(r, C) and n(r, £) is the origin of the mode mixing phenomena in |46j ]. 

The inclusion of a tensor matter component on the brane will source massive modes. The main source of this will 
be the photon quadrupole which develops late in the history of the universe. This issue will even more important in 
the case of scalar perturbations since the perturbations in the density, pressure and velocity will contribute to this. 
Some of the methods discussed here should be of use in taking this effect into account. 
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APPENDIX: ORTHOGONALITY RELATIONS 



Sturm-Liouville theory tells us that our basis functions will be orthogonal, but we need to know the normalization. 
We will derive the orthogonality here explicitly for the Bessel and parabolic cylinder functions used and determine 
the coefficient of the delta function in the orthogonality relation. 
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There is a standard formula for Bessel functions which tells us that 

xD u (mx)D l ,(nx) dx = — - (mD^+i(mx)D L ,(nx) — nD^imx^Du+iinx)] , (A.l) 

where D v and D v are any two Bessel functions of order v and m =/= n. This can be derived by an argument given 

in jljj • ^ we consider two functions g± (x) and g 2 (x) which satisfy the differential equations <?" + P n (x)g n = for 
some functions P n , then 

/<*-*>»■**- w& ("> 

Applying this to Bcssel's equation, we see that 

xD v (rnx)D v (nx) Ax oc S(m — n) , (A. 3) 



if D v {nx) and D v (nx) are chosen to vanish, or have vanishing derivatives, on the boundary. We are interested in the 
specific case 

D v (mx) = C^m-jx) = Yi(ml)J v {mx) — Ji(ml)Y^(mx) , (A. 4) 

with v = 1, and 



Cu(m;x) ~ \ — — (Yi{ml) 2 + Ji(mZ) 2 ) sin (mi + phase angle) , (A. 5) 

V Trmx \ J 

asymptotically, when x is large. So the constant of proportionality can be deduced ^lj to be 

/ D ly (mx)D u (nx )dx = — (Yi(mZ) 2 + Ji[ml) 2 ) 5(m-n). (A.6) 
Jxq m\ J 

For the time-dependent parts of the mode solutions, ipm(T) — r~ 1 j/ m (r), we can construct a similar orthogonality 
relation The y m satisfy the differential equation 177|l 

d 2 y m i ( ^ i m 2 A 2 2 



dt 2 \ k A 
so applying the same argument as before gives 



Hh = 0, (A.7) 



A 2 f 

-j^(m 2 ~ n 2 ) I t 2 y m y n dt = y m y' n - y' m y n , (A. 



where the r.h.s. will be zero due to the boundary conditions. To evaluate J t 2 y m y n dt when m = n, we can use the 
asymptotic form H8UI) 

e k 2 -n/8mA f j ^2 ~| 

Vm ~ (W)VV7 ° XP {-2 mAr2 ~ W l0gT + l6m ) ' (AJ) 

to write the integrand as 

e k^/SmA e k 2 n/SnA f f if ,2 - fc 2 



■ kt exp |-|(m + n)Ai 2 - log^" 1 *) - i^L log(fc~4) J (A.10) 



V2A(mny/ A 

Neglecting the log terms in the exponential and making the change of variables x = At 2 



k 2 -it / 8mA k 2 -n / 8nA 

t 2 y m y n dt~ * ,„ rrji i e' lmx ' 2 e~ mx ' 2 dx (All) 



(2^1) 3 / 2 (mn) 1 /4 

which is asymptotically the same as J xJ )Ji (mx/2)J ll (nx/2), allowing us to calculate the coefficient of the delta function. 
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